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A  Survey  of  Knowr)  Indicators  of  Auroral  Substorm  Onset 


I.  INTRODUCTION 


The  study  of  when  and  where  auroral  phenomena  will  occur  has  become  an  area 
of  major  concern,  particularly  in  recent  years.  As  scientists  learn  more  about  the 
mechanisms  for  various  auroral  emissions  through  repeated  measurements,  the 
measurements  themselves  have  become  increasingly  more  refined.  As  the  interest 
in  the  more  subtle  emissions  increases,  scientists  are  seeking  to  measure  intense 
discrete  auroral  emissions  of  Class  III  aurorae  or  greater. 

Early  ground  observations  usually  involved  little  more  than  a  knowledge  of 
periodic  solar  activity  and  ground  based  magnetometer  variations.  Once  a  substorm 
commences  it  can  be  immediately  measured  from  the  ground.  The  evolvement  to 
more  complicated  airborne  measurements  requires  extensive  planning,  both  in  the 
iong  term,  due  to  increased  costs,  and  in  the  short  term,  due  to  delays  between 
launching  the  vehicle  and  acquiring  the  often  erratic  auroral  event. 

Short  term  predictions  are  useful  for  experiments  on  vehicles  capable  of  dynamic 
measurements  such  as  the  space  shuttle,  which  may  prove  to  be  one  of  the  most 
effective  vehicles  for  auroral  acquisition.  The  shuttle  has  the  capability  of  acces¬ 
sing  large  sections  of  the  auroral  zone  from  an  altitude  of  approximately  200-400  km. 
Given  indications  within  30  min  as  to  the  region  for  the  highest  likelihood  of  auroral 
occurrence,  the  orbiter  can  be  oriented  such  that  it  could  acquire  the  auroral  event. 
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The  approach  of  this  report  will  be  to  discuss  two  kinds  of  auroral  indicators. 
The  first  section  discusses  long  term  indicators.  These  involve  solar  events  such 
as  the  sunspot  cycle,  solar  flares  and  coronal  holes.  Observations  of  this  type 
usually  involve  lead  times  of  one  day  to  approximately  eleven  years. 

The  second  section  discusses  short  term  indicators,  with  lead  times  ranging 
from  instantaneous  to  an  hour  or  so  before  substorm  onset.  These  indicators 
involve  satellites  and  ground  measurements.  There  are  periodic  measurements 
of  the  auroral  oval  by  low  orbiting  (~800  km)  satellites  that  pass  over  roughly  the 
same  geomagnetic  zone  once  every  102  minutes.  They  show  the  location  of  the 
equatorward  boundary  of  the  oval  and  the  general  level  of  substorm  activity  in  the 
oval.  In  addition  there  are  higher  satellites  at  geosynchronous  orbits  (about  6.6 
earth  radii )  which  can  measure  particle  flux  densities  and  magnetic  fields,  and 
even  higher  satellites  which  can  measure  the  intensity  of  the  interplanetary  mag¬ 
netic  field  in  the  bowshock  region  outside  the  magnetosphere  towards  the  sun. 

(There  are  no  such  operational  satellites  in  this  region  in  the  present  year,  1985. ) 
Ground  based  measurements,  such  as  all  sky  visual  observations  and  magnetograms, 
are  useful  for  tracking  the  phase  of  the  substorm. 


2.  LONG  TERM  INDICATORS 


2.1  The  1 1-Year  Suiupol  Cycle  and  the  27-Day  Recurrence  Tendency 

All  considerations  for  the  source  of  auroral  electrons  and  protons  begin  with 

the  sun.  Observations  of  sunspot  activity  have  indicated  a  definite  periodicity  of 

eleven  years.  This  period  is  called  the  sunspot  cycle.  During  the  cycle,  sunspot 

activity  generally  rises  to  a  maximum  and  then  subsides.  It  is  between  the  maximum 

and  the  end  of  the  cycle  when  auroral  activity  is  known  to  be  highest.  According  to 

Chamberlain^  it  is  about  two  years  after  this  sunspot  maximum.  But  recent  studies 
2 

by  Gneryshev  have  indicated  that  the  maximum  is  actually  the  first  of  two  maxima, 
the  second  of  which  occurs  about  two  to  three  years  after  the  first.  Gneryshev 
characterizes  the  first  maxima  by  a  maximum  number  of  small  sunspots  and  flares, 
a  homogeneous  corona,  and  few  auroral  occurrences.  The  second  maxima  is  charac¬ 
terized  by  bigger  sunspots  and  flares,  a  very  inhomogeneous  corona,  and  a  maximum 
number  of  auroral  occurrences. 


1.  Chamberlain,  J.  W.  (1961)  Physics  of  the  Aurora  and  Airglow.  Academic  Press, 

New  York. 

2.  Gneryshev,  M.  N.  (1977)  Essential  features  of  the  11-year  solar  cycle. 

Solar  Physics,  51:175-183. 
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The  sun  rotates  with  respect  to  the  earth  once  every  27  days.  Sources  of  high 
speed  streams  in  the  solar  wind  on  the  sun's  surface  tend  to  rotate  as  rigid  volumes 
with  only  a  very  slight  variation  of  their  rotational  period  with  latitude.  ^  Hence 
there  is  a  tendency  for  the  effects  of  these  sources  to  repeat  every  27  days.  Musical 
diagrams  of  the  Kp  index  are  sometimes  suggestive  of  this  recurrence  tendency, 
depending  on  the  stage  of  the  sunspot  cycle.  One  such  diagram  is  shown  in  Figure  1 

A 

(after  Mayaud  )  for  the  year  1933  when  recurrence  was  high. 
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Sargent  has  combined  the  27  -day  recurrence  phenomenon  in  the  aa  index 
(instead  of  Kp)  with  the  overall  periodicity  of  the  sunspot  cycle  every  11  years,  the 
results  of  which  are  shown  in  Figure  2.  The  aa  index  is  based  upon  the  basic  K 
indices  and  is  probably  one  of  the  best  indicators  of  long  term  geomagnetic  activity 
based  on  the  homogeneity  of  the  series  and  the  indices'  reliability  as  a  worldwide 
indicator.  ®  In  his  analysis  Sargent^  has  computed  a  correlation  coefficient  by 
matching  a  set  of  54  aa  indices  for  a  27  -day  rotation  with  the  set  of  54  aa  indices 
for  the  next  27  -day  rotation  and  assigned  this  correlation  coefficient  for  the  whole 
54 -day  period. 


Figure  2.  Twenty -Seven-Day  Recurrence  Index  (solid  line)  and  Smoothed 
Monthly  Mean  Sunspot  Numbers  (dashed  line)  From  1968  to  Present. 
Sunspot  cycles  are  numbered  according  to  convention 


He  noted  the  following  features  in  the  trends  of  recurrent  activity; 

(1)  The  onset  of  the  recurrent  activity  is  usually  abrupt,  and  for  much 
of  the  past  116  years  it  is  possible  to  say  confidently  whether  or 
not  the  signal  is  present. 


5.  Sargent,  H.  H,  ( 1985)  Recurrent  geomagnetic  activity:  evidence  for  long-lived 

stability  in  solar  wind  structure,  J,  Geophys,  Res.  ,  ^;1425-1428. 

6.  Mayaud,  P.  N.  (1973)  A  hundred  year  series  of  geomagnetic  data  1868-1967, 

lAGA  Bull,  Vol.  33, 


(2)  The  termination  of  the  recurrent  period  is  almost  always  as  abrupt 
as  the  onset  and  occurs  about  1  year  into  the  new  sunspot  cycle. 

(3)  The  duration  of  the  recurrent  phases  is  longer  in  some  cycles 
than  in  others.  These  longer  phases  occur  every  22  years.  The 
lengths  of  these  periods  of  recurrent  behavior  are  surprisingly 

long,  accounting,  at  times,  for  as  much  as  5  years  of  a  sunspot  cycle. 

2.2  Coronai  Hole* 

Coronal  holes  are  regions  of  unusually  low  density  and  low  temperature  in  the 
solar  corona.  It  has  generally  been  accepted  that  the  27  -day  recurrent  behavior  is 
directly  related  to  these  open  field  structures  on  the  sun.  ^  which  are  sources  of 
high  speed  streams  in  the  solar  wind.  This  acceptance  has  been  based  upon  observa¬ 
tions  that  long  lived  patterns  of  coronal  holes,  high  speed  solar  wind  streams,  and 
geomagnetic  disturbances  appear  remarkably  alike.  Sheeley  has  noticed  that  the 
polarity  of  the  magnetic  field  emerging  from  the  base  of  a  coronal  hole  is  the  same 
as  the  predominant  polarity  of  the  interplanetary  magnetic  field  carried  by  the 

O 

corresponding  high  speed  stream. 

Sheeley  found  the  transit  time  from  passage  of  the  coronal  hole  over  the  sun's 
central  meridian  to  the  time  of  geomagnetic  disturbance  to  be  about  2.  5  days, 
assuming  the  latitude  of  the  hole  is  close  to  the  heliographic  latitude  of  the  earth. 
However,  this  is  only  an  average  figure.  Sometimes  the  streams  arrive  at  the  earth 
before  the  hole's  central  meridian  passage  or  after.  This  may  be  due  to  the  fact 
that  the  boundaries  of  coronal  holes  sometimes  extend  non -radially  from  their  bases 
in  the  chromosphere  into  the  outer  corona.  Koomen  also  deduced  from  satellite 
measurements  at  3-10  Re  that  photospheric  magnetic  flux  emerging  from  low  lati¬ 
tudinal  coronal  holes  may  fan  out  considerably  with  distance  from  the  sun.  Sheeley 
suggests  the  low  latitude  flux  could  join  with  flux  from  a  polar  coronal  hole  to  form 
a  large  scale  interplanetary  magnetic  sector.  He  concludes  that  "since  coronal 
holes  are  visible  from  the  earth  several  days  before  they  are  carried  across  the 
central  meridian  by  solar  rotation,  and  since  2-3  days  remain  before  the  wind  from 
the  hole  would  reach  the  earth,  we  should  expect  to  be  able  to  predict  the  arrival 
of  the  high  speed  streams  and  their  associated  magnetic  effects  approximately  a  week 
in  advance. " 


7.  Suess,  S.  T.  (1979)  Models  of  coronal  hole  flows.  Space  *  jience  Review,  ^^;lS9-2()0, 

8.  Sheeley,  N.  R. ,  Harvey,  J.W. ,  and  Feldman,  W.C.  ( 197 6)  Coronal  holes,  solar 

wind  streams,  and  recurrent  geomagnetic  disturbances:  1973-1976, 

Solar  Physics,  ^;27 1-278. 

9.  Koomen,  M.  J. ,  and  Howard,  R.  A.  (1974)  Solar  Physics,  j7:469. 


2..‘t  Solar  Klarvs 


A  solar  flux  is  a  sudden  brightening  of  the  photosphere  that  remains  visible 
lor  several  hours.  Solar  flares  tend  to  occur  near  sunspots  and  their  number  can 
b»'  related  to  the  number  of  sunspots*®  by 

N|.  =  <v(R-  10)  . 

Where  Nj.  is  the  number  of  flares  per  solar  rotation,  is  the  mean  sunspot 

tuimber  and  a  is  a  constant  value  between  1.5  and  2. 

The  most  spectacular  aurorae  are  usually  associated  with  solar  flares.  * 

Newton*  '  found  that  the  largest  flares,  when  they  occurred  within  45*  of  the  center 

if  the  disk,  wc-re  followed  in  80  percent  of  the  cases  by  a  magnetic  storm  within  a 

maximum  lime  of  three  days.  The  most  likely  time  interval  between  a  great  flare 

12 

and  its  associated  magnetic  storm  is  about  one  day  or  slightly  less,  but  many 

occur  one  to  two  days  after  the  event.  Generally,  stronger  flares  have  shorter 

lead  times  before  the  substorm  onset  than  the  lead  times  of  weaker  flares. 

13 

Dodson  and  Ilederman  '  showed  that  a  flare  was  especially  effective  (92  percent 
of  the  cases  studied)  if  it  is  accompanied  by  solar  radio  bursts  below  200  Mhz 
before  the  flare  reaches  maximum  brightness. 

Krivsky*'^  did  a  similar  study  but  only  analyzed  solar  radio  bursts  at  the  fre¬ 
quencies  of  536  and  606  MHz,  and  with  200  solar  flux  units  (solar  flux  unit  = 

•22  2  -1 

10  W  m  Hz  ).  Krivsky  first  compared  proton  flare  occurrence  with  Ae  index 
levels  for  32  superimposed  epochs.  The  patterns  showed  a  rise  in  the  Ae  index 
beginning  three  days  prior  to  the  flare  event  and  peaking  on  the  second  day  after 
the  flare,  followed  by  a  decrease.  A  similar  analysis  of  the  radio  bursts  at  538 
and  606  MHz  indicated  the  Ae  index  rose  3  to  4  days  prior  to  the  burst  and  peaked 
6  to  8  days  after  the  burst.  Krivsky  concluded  that  a  rise  in  the  Ae  index  could 
be  associated  with  an  active  sunspot  region  where  multiple  non -proton  radio  bursts 
occurring  before  and  after  the  flare  accounted  for  a  rise  in  the  average  Ae  index 
value  to  a  noise  level  of  about  230  gamma  prior  to  and  after  an  enhanced  peak  from 
a  major  flare  event.  These  results  are  shown  in  Figure  3, 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  39.  ) 
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Figure  3b.  Curve  of  the  Average  Values 
of  the  AE  Index  10  Days  Prior  to  and 
12  Days  After  the  Occurrence  of  the 
Solar  Radio  Burst  With  a  Flux  in  Excess 
of  200  s.  f.  u.  (solar  flux  unit  « 

-22  -2  -1 

10  W  m  Hz  on  a  Frequency  of 
536  MHz  (Ondrejov  Observatory). 
Number  of  epochs  152,  in  years 
1966-1973.  Left-hand  scale  is  in  units 
of  y.  The  curve  of  the  histogram  also 
gives  the  numbers  of  occurrences  of  the 
maxima  of  the  AE  index  on  the  individual 
days  (right-hand  scale),  the  data  being 
adopted  from  the  individual  epochs 
employed 


Figure  3c,  Curve  of  tlie  Average 
Values  of  the  AE  Index  10  Days  Prior 
and  12  Days  After  the  Occurrence  of 
the  Solar  Radio  Burst  With  a  Flux  in 
Excess  of  200  s.  f.u.  on  a  Frequency 
of  606  MHz  (Sagamore  Hill  Observatory) 
Number  of  epochs;  120  from  interval 
of  years  1966-1973.  Left-hand  scale 
is  in  units  of  y.  The  curve  of  the 
histogram  gives  the  number  of 
occurrences  of  the  maxima  of  the 
AE  index  on  the  individual  days 
(right-hand  scale),  the  data  being 
adopted  from  the  individual  epochs 
employed 


3.  SHORT  TERM  INDICATORS 


3.1  Recunrcnce  Probability  of  Auroral  Activity  in  DMSP  Imagea 

Much  has  been  learned  about  the  auroral  oval  from  the  DMSP  satellites,  which 

scan  a  path  approximately  3000  km  wide  on  the  earth's  surface.  The  altitude  of 

the  satellites  is  about  800  km.  The  mapping  takes  into  account  a  very  large  portion 

of  the  auroral  oval  so  that  a  fair  evaluation  can  be  given  to  the  level  of  auroral 
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activity  and  the  latitudinal  extent  of  the  oval.  Sheehan  et  al  rated  749  southern 

hemispheric  DMSP  images  taken  during  1972  and  1973  and  rated  them  according 

to  active  (A),  moderate  (M),  quiet  (Q)  and  no  aurora  (N>,  finding  the  overall 

occurrence  rates  of  each  to  be  0.  17,  0.55,  0.  16,  and  0.  1 1  respectively  for  those 

years.  He  also  found  that  after  one  DMSP  orbit  both  A  and  Q  tended  to  recur  in 

one -half  of  the  cases.  As  a  result  there  was  some  doubt  cast  on  Akasofu's  early 
16 

substorm  model  since  Sheehan  concluded  that  it  was  possible  for  magnetic 

17  18 

bays  of  several  hours  to  occur.  In  Figure  4,  Sheehan  geographically  shows  an 
example  of  how  the  entire  oval  substorm  should  be  considered  in  the  determination 
of  magnetic  bay  deviations  rather  than  by  the  duration  of  a  bay  at  a  particular 
station,  stating  "The  magnetic  signature  at  a  fixed  location  on  earth  develops  as 

Q 

the  earth  rotates  under  an  essentially  stationary  auroral  pattern. " 

18 

Sheehan  has  tabulated  the  number  of  occurrences  of  the  A,  M.  Q.  N,  and  X 
(missing  data)  ratings  with  respect  to  the  rating  from  the  previous  pass  and  has 
constructed  a  matrix  called  a  persistence  table. 

Two  such  tables  are  shown  in  Table  1  and  Table  2.  Sheehan's  persistence 
tables  indicated  significant  likelihood  for  activity  levels  to  recur  on  the  next 
pass  102  min  later. 


15.  Sheehan,  R.  E. .  Carovillano,  R.  L. ,  and  Gussenhoven.  M.  S.  (1982) 

Occurrence  and  recurrence  in  auroral  activity  in  DMSP  images, 

J.  Geophys.  Res. .  ^(No.  A3 ):358 1-3589. 

16.  Akasofu,  S.  -I.  (1968)  Polar  and  Magnetos pheric  Substorms,  D.  Reidel, 

Hingham,  Massachusetts, 

17.  Carovillano,  R.  L. .  Siscoe,  G.  L. .  Sheehan,  R.  E. ,  Eather,  R.  H. .  and 

Gussenhoven,  M.  S.  (1975)  Unified  Model  of  Auroral  Substorm 
Development.  AFCRL-TR-75-055»,  AUAdzbsIs. 

18.  Sheehan.  R.  E.  Ph.  D.  Thesis  (Boston  College,  1977),  p.  182. 


Figure  4.  Schematic  Representation  of  the  Relationship  of  Magnetic  Bay 
Activity  and  Overhead  Aurora.  Xq  and  denote  lines  of  constant  geo¬ 
graphic  and  geomagnetic  latitude,  respectively.  The  shaded  zones  denote 
areas  of  active  aurora  in  a  DMSP  image.  The  upper  inset  gives  the 
horizontal  (N-S)  component  of  the  magnetic  field  at  Station  1  as  it  rotates 
in  UT  under  the  aurora.  The  lower  inset  traces  the  horizontal  component 
of  the  magnetic  field  along  at  the  fixed  UT  which  applies  to  the  DMSP 
image  (after  Sheehan,  1977) 


In  addition  Sheehan  has  developed  a  scheme  for  computing  the  probability  that 
an  image  rating  will  be  followed  by  the  same  or  another  image  rating  for  any  number 
of  orbits  later.  Defining  Pj^  to  be  the  occurrence  frequency  of  an  image  rated  at 
index  j  to  be  immediately  followed  by  one  oflndex  i,  the  Joint  probability  that  an '  ' 
image  of  index  j  is  followed  by  one  of  index  i  two  orbits  later  is  given  by 


Table  1.  Persistence  Table  for  August  1972.  Elements  in  table  give  the  number 
of  cases  and  percentage  (in  parentheses)  that  images  with  a  chosen  rating 
(column)  were  followed  by  images  with  the  same  or  any  other  rating  one  orbit  later 
(row).  For  example,  the  element  9(18)  in  column  A,  row  M  indicates  that 
9  images  (18  %)  of  the  49  rated  A  in  August  1^72  were  followed  by  M  rated 
images  one  orbit  later  (after  Sheehan.  1977)^" 


A 

M+ 

August  1972 

M 

M- 

Q 

N 

A 

20(41) 

9(27) 

7(15) 

5(11) 

2(04) 

2(05) 

M+ 

7(14) 

5(15) 

8(17) 

3(07) 

3(05) 

0(00) 

M 

9(18) 

9(27) 

8(17) 

5(11) 

7(13) 

3(07) 

M- 

5(10) 

3(09) 

7(15) 

9(20) 

7(13) 

3(07) 

Q 

3(06) 

3(09) 

7(15) 

12(27) 

9(16) 

9(21) 

N 

0(00) 

2(06) 

1(02) 

3(07) 

8(15) 

16(37) 

X 

5(10) 

2(06) 

9(19) 

7(16) 

19(35) 

10(23) 

Total 

49(100) 

33(100) 

47(100) 

44(100) 

55(100) 

43(100) 

Table  2.  Combined  Persistence  Table  Based  on  All  Images  (748)  in  August  1972 
and  July  -  August  1973  (after  Sheehan.  1977)^® 


A 

M+ 

Combined 

M 

M- 

Q 

N 

A 

52(50) 

20(24) 

26(14) 

5(04) 

4(03) 

2(02) 

Mf 

23(18) 

9(11) 

21(11) 

3(02) 

3(02) 

0(00) 

M 

19(15) 

26(31) 

61(32) 

23(17) 

14(11) 

4(05) 

M- 

6(05) 

6(07) 

22(11) 

32(24) 

21(17) 

8(10) 

Q 

4(03) 

3(04) 

19(10) 

30(22) 

19(15) 

13(18) 

N 

0(00) 

2(02) 

6(03) 

9(07) 

15(12) 

25(30) 

X 

27(21) 

17(20) 

37(19) 

34<2») 

47(38) 

3K37) 

ToUl 

131(100) 

83(100) 

192(100) 

138(100) 

123(100) 

83(100) 

The  indexing  process  is  diagrammatically  shown  in  Figure  5. 
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Figures.  Sheehan,  Ph.  1).  Thesis,  (1977)18 


The  procedure  can  be  generalized  for  any  number  of  orbits  (n) 
n-1 


ik 


Ij  Jk 


If  the  process  is  extended  to  enough  orbits  the  persistence  characteristics 
would  eventually,  at  least  for  the  data  set  for  1972-1973,  return  to  the  overall 
recurrence  frequencies:  namely,  A  »  0.  17,  M  =  0.  55,  Q  «  0.  16  and  N  ■  0.  11 

For  experimental  purposes  the  initiation  of  the  construction  of  a  persistence 
table  perhaps  a  month  or  more  in  advance  of  a  high  altitude  auroral  measurement 
experiment  could  be  a  useful  tool  in  the  evaluation  of  the  persistence  of  substorm 
activity,  since  the  overall  activity  would  be  dependent  on  recent  solar  activity. 


3.2  The  Solar  Wind  and  the  Interplanetary  Magnetic  Field 

19 

It  was  first  suggested  by  Oungey  that  a  southward  turning  of  the  interplanetary 
magnetic  field  (IMF)  between  the  sun  and  the  earth  might  affect  levels  of  geomag¬ 
netic  activity  and  various  solar  wind  parameters  such  as  the  wind  velocity,  the 
southward  component  of  the  IMF  (B^),  and  the  total  strength  of  the  IMF  (B). 

Table  3  shows  some  analytical  expressions,  consisting  of  these  parameters,  which 

20 

correlate  fairly  well  with  geomagnetic  activity.  Bather  however,  has  warned 


19.  Dungey,  J.  W.  (1961)  Interplanetary  magnetic  field  and  the  auroral  sones, 

Phys.  Rev.  Letters,  ^47. 

20.  Bather,  R.  H. ,  Mende,  S.  B. ,  and  Weber.  E.J.  (1979)  Dayslde  aurora  and 

relevance  to  substorm  current  systems  and  dayside  merging, 

J.  Geophys.  Res. ,  84:3339. 


that  direct  correlation  of  geomagnetic  activity  with  can  be  misleading;  he 

suggests  that  B  changes  are  not  the  causative  driving  mechanism  but  that  sustained 

^  21 
|)eriods  of  negative  B^  (S  1  hour)  increased  the  probability  of  substorms.  Eather 

has  found  better  correlation  of  Ae  with  the  dayside  cusp  equatorward  boundary 
(directly  associated  with  substorm  activity)  than  with  B^  (see  Figure  6).  The  ac¬ 
curacy  of  the  Ae  index  as  a  measure  of  substorm  activity  is  about  0.7  according 
22 

to  Kamide  (see  Figure  6).  The  time  of  increased  probable  onset  of  a  substorm 

from  the  time  the  solar  wind  impacts  the  magnetosphere  has  been  estimated  to  be 

22  23  24 

about  one  hour  by  Kamide.  although  Baker  et  al  and  Meng  thought  it  to  be 
more  like  40  minutes. 


Table  3.  List  of  the  Correlation  Studies  Between 
Geomagnetic  Indices  and  Solar  Wind  Parameters 


Geomagnetic  iintet 

Rclaiiott  with  solar 

(time  resolution) 

wind  parameters 

Sn>d<r  w  jI  ( l**6.M** 

ZKp 

(Nhr) 

lilfp  -  {V  -  )30)/S.44 

OIbcrt  aa 

ZKp 

(7*  hr) 

ZKp  *  t>  -  262)/6.3 

B.llir.f  j/  (1464)  ” 

Kp 

(3  hr) 

..,(-5^)1 

Arnoldy  (1471)” 

AS 

(Ihr) 

AE  -  -O  lHEB.rt^  -  0.9l(U,t). 

•  -0)3(rp,r)g4>  0.12^0 

Bobrov  (H73)  *• 

Kp 

(3  hr) 

«)>./().  »,.j»,) 

Gjrrcu  w  o/  (N74)**  31  up,  AS 

(Ihr) 

ap.  AE  <  c, .  Cjrd.)  ♦  c.r* 

Murjtdtnd  dfitl  Hakimada  (H75| 
Burton  w  a/.  (IVTj)  M 

AE 

(Ihr) 

AS  •  Cf.K* 
d 

Du 

(2.5  mm) 

--Diu  •  FiSi^aDu^ 
di 

Crook»r  ti  a/.  (1977)** 

.ApaiSo  IO-*i,F*-l.9 

Sv.l»jafd(l977)  ** 

am 

(3  hr) 

1.17 

OaSm***)"* 

Maciawa  (1979) 

Ai.  AV.  am 

(Ihr) 

M 

Holier  and  Slavin  (1979) 

AL 

(Ihr) 

Rale  of  croiMm 

Raitofrftani 

Mar^aina(l979)  ^ 

AL 

(Ihr)  tt-«|*.*e.5)K»ni.S,)p(*‘’ 

Hura)afna«r«f.  (I9BI  i* 

AL 

(1  hr)  ALiy*  >  X(S,  »  as) 

Bak««ai:(l9tl)>  «« 

AS 

(2.5  mill)  rs,.r>s.  Sir,  MC. 

Ui<f««.((9lir** 

C«mlaiioii  ht«»M6  Um  hiuipliiimtii  mt  iM|MieiplHt(c  ctocuic 

IW* 

atiSMa(.(l9tl)*  ** 

Comlaiion  hatofan  * 

Hotter  and  Stavia(l9l2r 

AL 

sr‘.s,r,s;r 

lijima  and  foumra  (1912)*  ^ 

Cdrralaiion  heiwcan  BaW^atifnad  carram  damiip  and  Rf,  6  f ’5^. 

etc. 

McloiiiM«(.(l9l2r  ** 

AE 

.Volf:  Fof  of  the  wprtitipin  (kUniiioMof  wch  wotniaw.  ttt  liw  ftOwtnct  in  Afcawfc  ^ 

The  papers  with  *  ere  refcitneed  in  this  paper. 


(Duelo  the  large  number  of  references  cited  above  in  the  text  and  in  Table  3,  they 
will  not  be  listed  here.  See  References,  page  39. ) 
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Perreault  and  Akasofu  have  recently  developed  the  parameter  £  as  a  measure 
of  the  power  generated  by  the  solar  wind -magnetosphere  dynamo,  drawn  from 
Akasofu's  theory  that  the  kinetic  energy  eventually  reaching  the  polar  ionosphere  is 
driven  by  the  solar  wind  —  IMF  disposition  and  is  not  related  to  unloading  processes 

internal  to  the  magnetosphere.  (This  view  is  disputed  by  Schindler^®  and  Heikilla. 
45 

Perreault  and  Akasofu  define  e  as 

,,„2  .40,2 
£  -  VB  sin  -K  1 
2  o 

where  V  is  the  solar  wind  velocity,  B  is  the  magnitude  of  the  IMF,  1^  denotes  the 
boundary  of  the  magnetopause  on  the  dayside  (taken  by  Akasofu  as  fixed  at  7  Re), 
and  0 


tan‘^  <By/B^) 


B  >  0 
z 


*  180  -  tan’  (B  /B  )  B,  <  0  . 

y  z  2 

23  43 

Baker  et  al  and  Meloni  et  al  performed  cross  correlation  analysis  of  the 
various  solar  wind  quantities,  including  £,  with  the  Ae  index.  Baker's  results  are 
shown  in  Figure  7  and  give  information  on  the  lead  times  between  surges  of  south¬ 
ward  Bjj  and  resulting  Ae  behavior.  Meloni's  correlations  of  solar  wind  parameters 
vs  Ae  are  presented  in  Table  4.  Meloni's  correlation  differ  from  Baker's  in  the 
fact  that  they  consider  correlations  of  Ae  with  £  for  all  B^  (southward  and  northward) 

as  well  as  correlations  of  Ae  with  £  for  southward  B  . 

z 


45.  Perreault,  P. ,  and  Akasofu,  S. -I.  (1978)  A  study  of  geomagnetic  storms, 

Geophys.  J.  R,  Astron,  Soc. ,  ^;547. 

46.  Schindler,  K.  (1976)  Solar  Physics,  47;91. 


47.  Heikkila,  W.  F.  (1983)  The  reason  for  magnetospheric  substorms  and  solar 

flares.  Solar  Physics,  88:329-336. 
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Figure  7a.  Correlations  of  the  AE  Index  With 
Solar  Wind  Velocity  Shown  Independently  for 
Data  Acquired  in  Individual  Months  From 
November  1973  Through  June  1974.  The 
ordinate  is  the  correlation  coefficient  and  the 
abscissa  is  the  number  of  minutes  by  which 
the  AE  samples  lagged  the  solar  wind  velocity 
samples 
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Figure  7b.  Correlations  of  the  AE  Index  With  the 
GSM  z -Component  of  the  Interplanetary  Magnetic 
Field 
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Table  4.  Summary  of  the  Cross  Correlation  Coefficients  Obtained 
in  the  Five  Magnetic  Storm  Intervals  for  Contemporary  Values  of 
AE  and  Interplanetary  Parameters  and  for  the  1-h  Delay  in  AE,  for 
Data  During  All  Values  and  for  Data  During  Only  B  South  Values 


3 -to  NOV  73 

22-29  JAN' 74 

16-23  APN‘74 

7-!4  JUN'  74 

3-10  JUL'7« 

■ 

AE  AE(«1) 

AE  AEU1) 

AE  AE(-1) 

At  AE(>1) 

AE  AE(*I) 

log  e 

0  3Z 

0  29' 

0  52 

0  60 

0  51 

0  30 

0  37 

0  47 

0  47 

032 

1 

SIN*  y 

0  44 

0  43 

0  32 

0  35 

0  28 

0  34 

0  60 

0  68 

063 

0  73 

-0  38 

-0  47 

_ 

-0  36 

-0  45 

-0  33 

-0  43 

-0  38 

-0  70 

•0  66 

-074 

mm 

-0  40 

-0  40 

-047 

-0  43 

-031 

-061 

-072 

-0  63 

-070 

log  c 

0  S8 

0  64 

0  72 

0  70 

0  66 

0  68 

061 

0  69 

0  62 

0  63 

1 

SIN«  ^ 

0  43 

0  37 

0  07 

0  03 

0  16 

018 

0  32 

0  39 

0  49 

0  33 

Bi 

-0  38 

-0  32 

-0  25 

-0  33 

-0  27 

-0  30 

-0  49 

-0  67 

-0  64 

-068 

V-Bi 

-0  49 

-0  60 

-044 

-0  49 

-0  49 

-0  49 

-039 

-0  78 

-0  60 

-0  62 

< 

0  38 

0  63 

_ 

0  53 

033 

0  56 

0  33 

0  33 

0  63 

0  35 

0  28 

3.3  The  Equatonvard  Bmuidary  of  the  Auroral  Oval 

3.  3.  1  OBSERVATIONS  OF  THE  RELATIONSHIP  OF  THE  EQUATORWARD 
BOUNDARY  TO  THE  DISCRETE  AURORAL  REGIONS 

The  auroral  oval  can  be  thought  of  as  the  projection  of  the  plasma  sheet  onto 
the  polar  ionosphere.  Figure  8  shows  geometry  of  the  earth  and  plasma  sheet 
within  the  magnetosphere.  The  electrons  and  protons  originating  in  the  magnetosphe 
tail  (below  and  to  the  right  in  the  figure)  are  driven  by  the  cross  tail  electric  field 
toward  auroral  latitudes  along  the  magnetic  field  lines  in  the  plasma  sheet. 

Satellite  observations  in  the  early  70's  indicated  that  there  were  two  types  of 

48 

auroras  that  occurred:  diffuse  and  discrete.  Aboard  the  ISIS-2  satellite  Lui  et  al 
analyzed  simultaneous  measurements  from  a  soft  particle  spectrometer,  energetic 


48.  Lui,  A.T.  Y. ,  Venkatesan,  D.  ,  Anger,  C.  D. ,  Akasofu,  S. -I. ,  Heikkila,  W,  J. , 
Winningham,  J.  D.  ,  and  Burrows,  .1.  R.  (1977)  Simultaneous  observations  of 
particle  precipitations  and  auroral  emissions  by  the  ISIS -2  Satellite  in  the 
19-24  MLT  Sector,  J,  Geophys,  Res. ,  82:2210. 


particle  detectors  and  an  auroral  scanning  photometer.  The  data  showed  that  the 

49 

diffuse  auroras  (<  1  Kr)  tended  to  occur  in  what  Winningham  referred  to  as  the 
central  plasma  sheet  and  the  discrete  auroras  corresponded  to  what  was  referred 
to  as  the  boundary  plasma  sheet. 


Figure  8.  Sketch  of  Magnetic -Field  Lines 
and  Various  Plasma  Regions  of  the  Mag¬ 
netosphere.  Only  the  northern  dusk 
quadrant  is  shown  for  simplicity.  Also 
drawn  are  contours  denoting  the  flow  or 
convection  of  plasma  in  the  equatorial 
plane.  All  magnetic -field  lines  shown  in 
this  figure  except  the  two  near  the  dusk 
meridian  lie  completely  in  the  noon -midnight 
meridian  plane.  The  Van  Allen  radiation 
belts  lie  within  and  just  beyond  the 
plasm  as  phe  re 


The  discrete  auroras  were  found  to  be  em  bedded  within  but  often  in  the  pole- 
ward  half  of  the  diffuse  aurora.  Figure  9  shows  the  relationship  for  quiet  and 
substorm  times. 


49,  Winningham,  J.  D,  ,  Yasuhara,  F. ,  Akasofu,  S. -I. .  and  Heikkila,  W.  J.  (1975) 
The  latitudinal  morphology  of  10  eV  to  10  KeV  electron  fluxes  during 
magnetically  quiet  and  disturbed  times  in  the  2100-0300  MLT  sector, 

J.  Geophys.  Res.,  80;3148. 
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Figure  9.  Schematic  Diagram  to  Illustrate 
the  Spatial  Relationship  Between  the  Plasma 
Sheet  and  the  Two  Types  of  Aurora. 
Characteristics  of  electron  differential 
energy  spectra  in  both  types  of  aurora  for 
quiet  and  substorm  times  are  also  illus¬ 
trated  (Lui  et  aP°) 


Lui  et  al^®  looked  at  three  kinds  of  auroras:  those  with  diffuse  auroras,  those 
with  bright  arcs  but  no  substorm  and  those  with  substorms,  noting  the  boundaries 
of  the  diffuse  and  the  discrete  regions  as  well  as  Kp.  These  results  are  shown  in 
Figure  10  and  Table  5.  This  data  suggests  that  if  the  equatorward  boundary  of 
the  diffuse  aurora  is  known,  along  with  the  amount  of  magnetic  activity,  the  loca¬ 
tion  of  the  region  of  discrete  forms  could  be  estimated, 

A  similar  analysis  of  the  relationship  of  the  discrete  and  diffuse  equatorward 

50 

boundaries  had  also  been  measured  by  Nagata  et  al  in  the  southern  hemisphere. 
These  results  are  shown  in  Figure  11.^® 


50.  Nagata,  T. ,  Kirasawa,  T. ,  and  Ayakawa,  M.  ( 1975)  Discrete  and  diffuse 
auroral  belts  in  Antarctica,  in  Report  of  Ionosphere  and  Space  Research 
in  Japan,  29:149-152, 


WM 


INVAAIANT  unruoc 
n  ro  « 


90  T%  TO  Vi 

^xUJL^ 

i*  •'•§»»•  *»%*'*>*«'^ 

•— •  •#*■*<*  9»  »f*«»'»#'|  4*‘-t9tm  t(«C*r«*  •»♦«♦•♦••••• 
cru-si  •••4  )*  4v*9*«« 

— ■*♦<10#'  J# 

— '■—  a«t<4«  0*  •»««■  •'«<♦'*•*•«• 

93. •'^f  3* 


I  l«  J-|‘,*J 


ii4 


■  .V'  f  ^  '  "■ 

I'fcvr  ^  :v>  ■' J 
sO.-*  1 


f  9>m9  HMWWff  ••  iiiitri 


Figure  10.  Schematic  Diagram  to 
Summarize  the  Results  of  Comparison 
Between  Regions  of  Auroras  and 
Particle  Precipitation.  Each  pass  is 
represented  by  three  lines  to  indicate 
the  extent  and  characteristics  of  par¬ 
ticle  precipitation  and  the  types  of 
aurora.  Precipicatlon  and  auroral 
emissions  in  the  polar  cap  are  not 
included.  The  numbers  on  the  right 
side  identify  the  events  as  given  in 
Table  5 


Table  5.  Passes  Selected  for  the  Comparison  of 
Measurements  on  Particle  Precipitation  and  Auroral 
Intensity  by  ICxperiments  on  Board  the  Isis  2  Satellite 


No. 

Date 

Tlae  Interval 
of  Paasea 

*.* 

MLT  of 
Projected 
Path 

Cate|ory  of 
Auroral 
Activity 

1 

Dec. 

11. 

1971 

0S20-0534 

Oo 

22.1 

B 

2 

Dec. 

12. 

1971 

0404-0417 

lo 

23.1 

B 

3 

Dec. 

13. 

1971 

0248-0303 

3^ 

23.4 

A 

4 

Dec . 

13. 

1971 

0442-0457 

2* 

22.7 

A 

5 

Dec. 

1^. 

1971 

0519-0533 

1- 

21.6 

C 

6 

Dec. 

15. 

1971 

0209-0221 

2o 

23.5 

A 

7 

Dec . 

16. 

1971 

0247-0301 

(K 

23.3 

A 

8 

Dec. 

16. 

1971 

0440-0456 

0+ 

21.4 

A 

9 

Dec. 

20. 

1971 

0325-0339 

Oo 

22.4 

C 

10 

Dec. 

20. 

1971 

0518-0532 

Oo 

21.8 

C 

11 

Dec. 

21. 

1971 

0208-0222 

!♦ 

22.9 

c 

12 

Dec. 

22, 

1971 

0440-0456 

34 

21.1 

A 

13 

Dec. 

26. 

1971 

0323-0340 

4- 

21.6 

A 

14 

Jan. 

4, 

1972 

0322-0339 

3o 

20.6 

A 

15 

Jan. 

10, 

1972 

0321-0336 

lo 

19.7 

B 

16 

Jan. 

11. 

1972 

0359-0415 

4e 

19.0 

B 

17 

Jan. 

11. 

1972 

0551-0605 

44 

18.5 

B 

*Th»  Kp  arc  taken  froa  the  3-liour  Interval  contalnlnf  the 
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Figure  11.  Average  Locations  of 
the  Poleward  and  Equatorward 
Boundaries  of  the  Discrete  Auroral 
Belt  (a)  and  the  Diffuse  One  (b)  as 
Dependent  on  Geomagnetic  ^  Index 
and  the  Geomagnetic  Local  Time 


Spiro  et  al^^  examined  data  from  the  Atmospheric  Explorer  C  and  D  satellites 
and  concluded  that  regions  of  the  most  intense  precipitation  corresponded  to  regions 
where  upward  field  aligned  currents  are  typically  observed  by  satellite  magnetometers. 


51.  Spiro,  R.W. ,  Reiff,  P.  H. .  and  Maher,  L.J, ,  Jr.  ( 1982 )  Precipitating  electron 
energy  flux  and  auroral  zone  conductances  -  an  empirical  model, 

J.  Geophys.  Res.,  87:8215-8227. 
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Figure  12  shows  the  current  flow  at  auroral  altitudes.  The  currents  are  typically 
found  near  the  poleward  edge  of  the  evening  auroral  zone  for  active  geomagnetic 
conditions  (Ae  >  300  y).  With  respect  to  local  time  distributions.  Spiro  concluded 
that:  (1)  the  region  of  most  intense  precipitation  for  a  given  magnetic  local  time  is 
located  at  higher  latitudes  in  the  premidnight  sector  than  in  the  post  midnight  sector, 
and  (2)  for  higher  geomagnetic  activity  there  is  a  tendency  for  the  evening  high 
latitude  flux  to  extend  into  the  early  morning  sector  and  the  morning  low  latitude 
peak  to  extend  westward  into  the  late  evening  sector.  Spiro's  contours  of  precipitating 
energy  flux  are  shown  in  Figure  13. 


2-18  kcv 
ELCCTSONS 


DISTANCE  PERPENDICULAR 
TO  ARC  (Km) 


Figure  12.  Model  Current  System  for  Auroral  Arc 
as  Fitted  to  Rocket  Magnetometer  Observations  by 
Park  and  Cloutier  (197  1).  The  vertical  currents 
of  Figure  3.  16A  are  ass\imed  to  be  so  remote  as  to 
produce  a  negligible  effect  (Courtesy  of  Journal  of 
Geophysical  Research;  Copyright  by  American 
Geophysical  Union) 
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Figure  13.  Contours  of  Constant 
Precipitating  Electron  Energy  Flux 
Sorted  According  to  Auroral 
Electro  jet  (AE)  Index.  The  outermost 
contours  correspond  to  a  bin -average 
-1  -2-1 

energy  flux  of  2.  5  X  10  erg  cm  s 
with  adjacent  contours  representing 
factor  of  two  increases  in  energy 
flux 


Pinning  down  the  specific  location  of  the  discrete  events  is  difficult  because  of 

the  rapid  changes  that  can  take  place  between  successive  satellite  passes.  Two  or 

three  satellites  passing  over  the  2100-0300  sector  can  help  considerably  to  monitor 

52 

the  locations  of  the  discrete  regions  as  the  substorm  progresses.  Rostoker  et  al 
describe  the  definition  of  this  aspect  of  the  substorm  as  follows: 

During  this  interval  there  may  be  a  sequence 
of  intensifications  of  the  westward  electrojet,  each 
associated  with  a  Pi  2  micropulsation  burst  and  a 
westward  travelling  surge.  As  the  substorm 
develops,  the  region  of  discrete  auroras  in  the 
midnight  sector  expands  poleward  and  westward 

(the  poleward  bulge) .  The  authors  on  the 

basis  of  their  experience,  feel  that  the  magnetospheric 
substorm  as  defined  above  has  a  time  scale  of 

4 

1-3  hours. 

18 

Sheehan's  "fixed"  auroral  pattern  (Sheehan  ),  suggests  that  a  high  inclination 
shuttle  orbital  track,  which  would  be  fixed  with  respect  to  the  sun  and  passing  over 
the  midnight  sector,  could  be  able  to  monitor  the  progress  of  the  substorm,  con¬ 
sidering  the  mirroring  effect  between  the  northern  and  southern  ovals. 


52.  Rostoker,  G. ,  Akasofu,  S.  -I. ,  Foster,  J. ,  Greenwald,  R.  A. ,  Kamide,  Y. , 
Kawasaki.  K. ,  Lui,  A.  T.  Y. ,  McPherron,  R,  L, ,  and  Russell,  C.  T. 
(1980)  Magnetospheric  substorms  -  definition  and  signatures. 

J.  Geophys.  Res.,  82(A4):1664. 

-  lAVV 
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3.  3.  2  MODELS  CONCERNING  THE  LOCATION  OF  THE 
EQUATORWARD  BOUNDARY 

One  of  the  earliest  attempts  to  classify  the  equatorward  boundary  was  by 
Feldstein  and  Starkov.  Using  a  broad  network  of  all  sky  cameras  from  61.  6  to 
75.  3  degrees  latitude,  the  oval  was  mapped  every  15  min  and  correlated  with 
the  Q  magnetic  index.  Their  determination  of  both  northern  and  southern  boundaries 
as  a  function  of  local  time  and  magnetic  activity  are  shown  in  Figure  14.  These 
determinations  have  come  to  be  known  as  the  "classical”  oval. 


Figure  14.  The  Dependence  of  the  Location  of  the  Northern 
and  Southern  Edges  of  the  Auroral  Belt  on  the  Indices  (a)  Q, 
(b)  Q  ,  and  (c)  K-  of  Magnetic  Activity  Around  Midnight 
(22-02  hours)  in  Local  Geomagnetic  Time.  0'  is  the 
corrected  geomagnetic  latitude.  The  dashed  line  represents 
the  center  line  of  the  auroral  belt 


54 

Eather  did  extensive  measurements  of  the  equatorward  boundary  of  the  day- 

side  aurora  through  the  use  of  a  keogram  camera  and  found  the  boundary  locations 

to  be  about  5^  higher  in  latitude  than  Feldstein  had  found  for  Q  ■  0.  The  keogram 

measurements  could  be  useful  in  characterizing  the  oval  in  general  on  a  real  time 

basis  by  relying  on  up  to  date  modeling  of  the  boundary.  Computer  codes  have  been 

developed  that  characterize  the  oval  based  on  the  input  of  a  geomagnetic  index 
55 

(Comfort  ).  Such  codes  can  be  easily  modified  to  accept  a  known  location  of  the 
boundary  for  a  particular  MLT  and  characterize  the  rest  of  the  oval  based  on  how 
this  location  correlates  with  previous  measurements  that  define  the  characterization. 


53.  Feldstein,  Y.  I,  ,  and  Starkov,  G.  V,  (1967)  Dynamics  of  auroral  belt  and  polar 

geomagnetic  disturbances.  Planet.  Space  Science,  15;209. 

54.  Eather,  R,  H.  ( 1984)  Oayside  auroral  dynamics,  J.  Geophys.  Res.. 

89:1695-1700. 

55.  Com)ort,  R.  H.  (1972)  Auroral  Oval  Kinematics  Program,  NASA  Report 

CR-61373. 
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5  6 

Gussenhoven  et  aF  ”  have  modeled  the  location  of  the  equatorward  boundary 
for  the  whole  oval  based  upon  DMSP  satellite  measurements.  The  model  accepts 
Kp  as  its  input  and  determines  the  boundary  from  the  equation 


r 


^  + 
o 


a  K 

P 


where  ^  is  the  boundary  predicted  from  the  regression.  Gussenhoven's  tables  of 
Q  and  a  are  shown  in  Table  6. 


Table  6.  /»  =  a  oK 
o  p 


MLT 

Number* 

^0 

2 

cct 

OUOO-OIOO 

312 

66.1 

-1.99 

-0.80 

OIOO  0200 

220 

65.1 

-155 

-0.68 

0400-0500 

267 

67.7 

-1.48 

-0.57 

050(VU600 

1123 

67.8 

-1.87 

-0.71 

0600-«)700 

2462 

682 

-1.90 

-0.74 

0700-0800 

3159 

68.9 

-1.91 

-0.76 

08(X)-O<MI0 

2159 

693 

-1.87 

-0.73 

(NOO  1000 

1178 

695 

-1  69 

-0.66 

1000-1100 

864 

69.5 

-141 

-0.57 

iioo-i:oo 

513 

70.1 

-1.25 

-0.52 

I’OO  l.too 

353 

69.4 

-0.84 

-0.3$ 

15(X>-1<100 

63 

709 

-0.81 

-0.34 

IbtIO  1700 

204 

71.6 

-1.28 

-0.66 

I7(X>-IX00 

526 

71.1 

-1.31 

-0.69 

WXPI'XIO 

997 

71.2 

-1.74 

-0.82 

I'XXi-’iXX) 

2469 

70.4 

-1.83 

-0.82 

:ooo:i(X) 

3309 

69.4 

-1.89 

-0.82 

2nx>-:^x) 

.3092 

686 

-1.86 

-0.79 

’:o(V2.KX) 

1482 

67.9 

-1.78 

-0.77 

2.MXX  24tX) 

461 

67.8 

-2.07 

-0.81 

‘Number  ol  boundaries  used  in  each  regression. 
tCorrelaiion  coedicient. 


Sauvaud  et  at  have  come  up  with  a  model  for  the 
the  22-24  MLT  sector  based  on  the  Ae  index  given  by 


equatorward  boundary  in 


56.  Gussenhoven,  M. S. .  Hardy.  D. A.,  and  Heinemann,  N.  (1983)  Syatematics  of 
the  Equatorward  diffuse  auroral  boundary,  J.  Geophys.  Res. , 

88:5692-5708. 


57.  Sauvaud,  J.  A. ,  Crasnier,  J. ,  Galperin.  Yu. I.,  and  Feldstein.  Y. I.  (1983) 

A  statistical  study  of  the  dynamics  of  the  Equatorward  boundary  in  the  pre¬ 
midnight  sector,  Geophys.  Res.  Letts. .  jy^749-752. 
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Inv.  lat  ((U-g)  =  A  ‘  H  In  (Se)  *  C  [ln{Se)| 


with  Ac-  iM-injt  the  avc-i-age  of  Ae  over  various  time  intervals.  The  constants  A,  B 

2 

anil  f  arc'  shown  in  Table  <  alonj;  with  the  square  of  the  correlation  (r  )  and  the 
stanilaril  (•rr-of  (Sli).  The  results  indicate  that  the  best  fits  occur  when  Ae  is 
aver-aticd  over  five-  hour  intervals  indicating  that  the  boundary  location  depends 
mainly  on  the  past  magnetic  activity. 


Table  7.  Hesult  of  the  Kitting  Procedure  for  Various  Average 
AK.  Values 


■raging  l-’oriod 

A 

B 

C 

2 

r 

SE 

20  m  inutes 

64.  81 

2.  61 

-0.  54 

0.77 

1.  67 

2  houfs 

60.  69 

4.  61 

-0.78 

0.  81 

1.  53 

3  lioufS 

59.  29 

5.  40 

-0.  88 

0.  86 

1.  33 

5  hours 

64.  11 

3.71 

-0.7  5 

0.  90 

1.  13 

10  hours 

57.  53 

6.41 

-1.02 

0.75 

1.74 

20  hours 

81.  11 

-3.  54 

0 

0.  61 

2.09 

3.4  Obaervaliom  of  Energetic  Electron  Intencity  Qiangee  at 
Geosynchronous  OrUt 


Studies  of  the  relationship  of  evening/midnight  sector  electron  dropouts  at 
geosynchronous  orbit  (about  6.  6  Re)  to  changes  in  the  auroral  electrojets  offer 
exciting  possibilities  for  the  prediction  of  the  onset  of  substorms  (Erickson  and 
Winckler^^,  Erickson  et  al^^,  and  Sauvaud  and  Winckler®®). 

Studies  of  the  behavior  of  energetic  electrons  (E  >50  KeV)  from  ATS-l,  ATS-5 
and  ATS-6  satellites  in  the  evening  and  midnight  sectors  have  shown  that  a  pro¬ 
nounced  decrease  of  energetic  electron  flux  regularly  precedes  by  about  an  hour 
the  marked  injection  to  synchronous  orbit  of  hot  plasma  from  the  tail  of  magnetosphere 

58.  Erickson,  K.  N. ,  and  Winckler,  J.  R.  (1973)  Auroral  electrojets  and  evening 

sector  electron  dropouts  at  synchronous  orbit.  J.  Geophys.  Res. .  ^:8373-8380. 

59.  Erickson,  K.  N. ,  Swanson,  R.  L, ,  Walker,  R.J. .  and  Winckler.  J.R.  (1979) 

A  study  of  magnetosphere  dynamics  during  auroral  electrojet  events  by 
observations  of  energetic  electron  intensity  changes  at  synchronous  orbit, 

J.  Geophys.  Res. ,  ^;93 1-942. 

60.  Sauvaud,  J.  A. ,  and  Winckler,  J.R.  (1980)  Dynamics  of  plasma,  energetic 

particles,  and  fields  near  synchronous  orbit  in  the  nighttime  sector  during 
magnetospheric  substorms.  Amer.  Geophys.  Union.  Report  No.  9A1S92. 
pp.  2043-2056.  -  - 


and  simultaneous  increase  of  the  flux  of  the  energetic  particles.  Coincident  with 

the  start  of  the  electron  decrease,  a  quiet  arc  system  appears  and  intensifies  as 

59 

the  decrease  and  field  inflation  continue  (Erickson  et  al  ).  During  some  decreases, 
activity  such  as  surges  or  brightenings  may  appear,  but  the  major  "breakup"  always 
coincides  closely  with  the  expansion  phase  onset.  Figure  15  is  a  representation  of 
this  process,  which  also  correlates  well  with  the  variations  of  the  electron  flux. 


Figure  15.  Schematic  Representation 
of  the  Simultaneous  Variations  of  the 
H  Component  of  the  Magnetic  Induction 
at  6.  6  R^  (antiparallel  to  the  dipole 
axis).  oAhe  AE  Index  and  of  the 
Energetic  Particle  Flux  at  6.  6  Rg 
During  a  Substorm  Preceded  by 
a  Quiet  Period 
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Variations  of  electron  fluxes  and  magnetic  field  variations  for  22  November 
1967  and  2  September  1987  aboard  ATS-1  are  shown  in  Figures  16  and  17,  along 
with  magnetometer  traces  for  a  number  of  ground  stations.  The  Au  and  Al  traces 
for  the  period  are  shown  as  well. 
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22  NOVEMBER  1967 

Figure  16.  Average  Electron  Counts 
Over  6-min  Intervals  in  Two  Energy 
Channels  for  22  November  1967.  The 
AU  and  Al.  indices,  magnetogram 
records,  and  time  scales  are  as 
described  for  Figure  2 
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2  SEPTEMBER  1967 


Figure  17.  Average  Electron  Counts 
Over  6-min  Intervals  in  Two  Energy 
Channels  for  2  September  1967, 

Lx>cal  B,  the  AU  and  AL  indices, 
magnetogram  records,  and  time 
scales  are  as  described  for  Figure  2 


In  the  evening  and  near -midnight  MLT  regions,  the  decreases  in  flux  of  the 
50-150  KeV.  150-500  KeV  and  500-1000  KeV  intervals  have  been  shown  to 
correlate  with  the  change  in  the  H  component  of  the  magnetic  B  field;  that  is.  when 
the  flux  decreases,  H  decreases,  and  when  the  flux  increases,  H  increases.  In 
the  morning  sector,  however,  some  large  variations  in  the  50-150  KeV  have  been 
measured  with  no  accompanied  change  in  the  B  field.  This  type  of  electron  flux 
increase  has  been  attributed  to  the  injection  of  plasma  sheet  particles  near  local 
midnight  followed  by  their  drift  toward  the  east  (Winckler;  *  Deforest  and  Mcllwain  ), 
where  they  arrive  somewhat  later  than  the  more  energetic  fluxes  150  KeV)  during 

the  expansive  phase  of  substorm.  They  follow  morphological  patterns  that  have 
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been  termed  as  "something  different  than  a  regular  substorm”  by  Mcllwain.  The 

location  of  the  injection  boundary  of  the  plasma  sheet  particles  at  substorm  onset 
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has  been  derived  by  Mcllwain  as 
Rjj  =  (122  -  10  Kp)  /  {<D  -  73) 

where  Rj^  is  the  injection  boundary  position  in  the  equatorial  plane  as  a  function  of 

K  and  MLT  in  <1800  S  ^  S  2400  CGLT). 

P 

Figure  18  shows  the  superpositioning  of  161  traces  of  the  Ae  index,  the  50-150  KeV 
fluxes,  the  150-500  KeV  fluxes,  and  the  B  field,  which  confirm  the  election  decrease 
during  the  growth  phase  of  the  substorm. 

The  geographic  location  of  a  substorm  can  be  roughly  approximated  with  mag¬ 
netic  field  line  models.  Erickson  has  successfully  monitored  the  growth  phase  of 
the  substorm  by  employing  these  field  line  models  in  combination  with  observations 

5 

of  the  changes  in  the  V  component  of  the  magnetic  field  aboard  the  GOES  satellites. 

At  the  present  time,  the  GOES  satellites  have  magnetometers  that  monitor  the  H 
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and  V  components  of  the  B  field  and  this  data  is  available  real-time  (Joselyn  >. 

61.  Winckler,  J.  R.  (1970)  The  origin  and  distribution  of  energetic  electrons  in  the 

Van  Allen  radiation  belt.  Particles  and  Fields  in  the  Magnetosphere, 

B.  M.  McCormac,  Ed. ,  D.  Reidel,  tlordrecht,  Netherlamds,  p.  332. 

62.  Deforest,  S.  E..  and  Mcllwain.  C.  E.  (1971)  Plasma  clouds  in  the  magnetosphere. 

J,  Geophys.  Res. ,  J[^3387 -361 1. 

63.  Mcllwain,  C.  E.  (1973)  Electron  beams  near  the  magnetic  equator.  Physics  of 

the  Hot  Plasma  in  the  Magnetosphere,  Plenum,  New  York. 


64.  Joselyn,  J,  A.,  and  Grubb,  R.  N,  (1985)  The  Space  Environment  Monitors 
Onboard  GOES,  Report  AIAA-83-0238.  AIAA  23rd  Aerospace  Sciences 
Meeting,  14-17  January  1985,  Reno.  Nevada, 


Figure  18.  One  Hundred  Sixty -One 
Superposed  Traces  of  AE  Index,  50- 
to  150 -keV  and  150-  to  500-keV 
Electron  Counts  and  Magnetic  Induction 
in  the  2100  2400  MLT  Bin.  Selected 
From  About  300  Days  of  Data  With 
KP  =  4+.  The  heavy  line  shows 
the  superposed  epoch  average  for  all 
161  events  in  this  local  time  bin. 

The  dots  with  dashed  lines  indicate  the 
background  curve  that  was  obtained  by 
a  superposition  of  all  the  300  24 -hour 
day  intervals  with  Kp  §  4+.  The  start 
of  each  event  (T  *  0)  was  chosen  at 
the  onset  time  of  the  150-  to  500 -keV 
electron  count  increase  [from  Swanson 
19781.  The  entire  24-hour  period  is 
shown  in  the  upper  panel  for  the 
same  T  =  0 


3.5  Ground  Bawd  Obwrvaliont 

3.  5.  1  VISUAL  ALL-SKY  CAMERA  AND  PHOTOMETERS 

All  sky  cameras  are  one  of  the  most  basic  means  to  spot  an  aurora  quickly 
while  carrying  on  other  ground  measurements  at  the  same  time.  Simultaneous  all 
sky  camera  pictures  from  various  ground  stations  can  begin  to  (riece  together  Uie 
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substorm  puzzle.  The  sequence  of  all  sky  pictures  shown  in  Figure  19  (Reference  65) 
from  several  meridian  chain  stations  show  a  considerably  better  picture  of  the 
substorm  phase  than  one  picture  does.  Unfortunately,  little  has  been  accomplished 
to  establish  a  central  facility  that  could  simultaneously  monitor  all  sky  camera 
pictures  from  all  ground  stations  on  a  real-time  basis,  although  this  has  been  done 
semi -real -time  from  aircraft  flying  along  the  oval. 


The  relative  amounts  of  the  constituents  identify  different  types  of  aurorae  which 
may  have  different  morphological  patterns. 

3.5.2  GROUND  BASED  MAGNETOMETER  TRACES 

Ground  based  magnetometer  traces  are  invaluable  for  identifying  a  geomag¬ 
netic  disturbance  in  the  vicinity  of  a  station  and  are  much  more  useful  than  the  Kp 
index  in  monitoring  geographically  the  phase  of  the  substorm.  Most  stations 
record  the  H.  D  and  Z  traces  which  are  the  components  of  the  total  magnetic  force 

field.  Figure  20  shows  the  geometry.  A  thorough  discussion  of  the  interpretation 
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of  the  H.  D  and  Z  component  traces  is  given  by  Rostoker  et  al.  Figure  21  shows 
a  general  example  of  how  these  traces  are  interpreted. 
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Figure  20.  The  Magnetic  Force,  F,  and  the 
Elements  of  the  Earth's  Magnetic  Field 
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Figure  21.  Schematic  Representation  of  Magnetograms  From  Various 
Sites  Placed  at  Different  Positions  With  Respect  to  the  High-Latitude 
Substorm  Current  Systems,  These  synthetic  records  can  be  compared 
with  real  data  to  establish  the  position  and  time  sequence  of  events  for 
an  interval  of  substorm  activity.  Typical  deflections  at  high-latitude 
observations  range  from  ~  100  to  1000  nT,  while  typical  deflections 
at  low -latitude  stations  range  from  a  few  nanoteslas  to  a  few  tens  of 
nanoteslas 


The  importance  of  numerous  stations  becomes  obvious  upon  examination  of 
one  component  trace  from  a  number  of  stations  along  a  meridian  chain.  An 
example  of  a  "stackplot"— the  superpositloning  of  the  traces  of  several  stations 
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top  of  one  another— is  shown  in  Figure  22  (see  Reference  66).  Such  chains  are 
able  to  monitor  the  expansion  and  contraction  of  the  oval,  and  pinpoint  the  intensity 
and  width  of  the  westward  or  eastward  electrojets. 
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Figure  22.  Stackplot  of  the  Z  Trace  From  the 
Alaskan  IMS  Magnetometer  Chain,  Illustrating 
a  Substorm  Signature  Near  0900  UT 


3.  5.  3  COORDINATED  GROUND -SATE  LUTE  AND 
AURORAL  RADAR  MEASUREMENTS 

An  important  phenomenon  to  be  considered  with  regard  to  the  onset  of  a  sub¬ 
storm  is  the  Harang  Discontinuity.  The  Harang  Discontinuity  is  a  thin  strip,  roughly 
aligned  from  east  to  west,  which  separates  the  eastward  and  westward  electrojets, 
and  its  behavior  is  strongly  dependent  on  magnetic  activity.  It  is  generally  recog¬ 
nized  that  the  first  brightening  arcs  of  a  substorm  onset  usually  occur  in  the  vicinity 
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of  the  Harang  Discontinuity  (Heppner;  Nielsen  and  Greenwald  ). 


66.  Joselyn,  J.  A.  et  al  (1982)  SESC  Geomagnetic  Prediction,  in  Proceedings  Work¬ 

shop  on  Satellite  Drag  (NOAA/ space  Environment  Laboratory, 

Boulder.  Colorado,  p.  135. 

67.  Heppner.  J.  P.  (1972)  The  Harang  discontinuity  in  auroral  belt  ionospheric 

.currents,  Geophys.  Publ. ,  ^9;  IPS. 

68.  Nielsen.  E. ,  and  Greenwald,  R.  A.  (1979)  Electron  flow  and  visual  aurora  at 

the  Harang  discontinuity,  J.  Geophys.  Res. ,  84:4189-4200. 


Nielsen  and  Greenwald  have  examined  the  Harang  Discontinuity  using  the 
STARE  (Scandinavian  Twin  Auroral  Radar  Experiment)  and  present  ionosphere 
electron  drift  patterns  which  were  associated  with  auroral  arcs  seen  by  all  sky 
cameras  (see  Figure  23).  They  noted  that  for  high  magnetic  activity  prior  to  the 
substorm  the  discontinuity  surged  equatorward  with  a  mean  speed  of  190  m/sec  to 
be  terminated  by  the  sudden  activation  of  an  arc  in  the  discontinuity,  followed  by  a 

poleward  expansion  at  substorm  onset.  The  arcs  tend  to  form  just  poleward  of  the 
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discontinuity  and  diffuse  aurorae  tend  to  occur  equatorward  of  it.  Heppner  and 
70 

Maynard  have  also  investigated  the  dynamics  of  the  discontinuity. 
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Figure  23.  The  Temporal  Evolution  of  the  Electron  Drift  Velocities  During 
High  Magnetic  Activity.  Note  the  equatorward  motion  of  the  Harang 
discontinuity.  In  this  case,  auroral  arcs  were  located  near  the  poleward 
and  equatorward  borders  of  the  discontinuity 


69.  Heppner,  J.  P.  (1976)  Aurora  and  Airglow.  B. M.  McCormac.  Ed., 

Reinhold,  New  York,  p.  75. 

70.  Maynard,  N.  C.  (1974)  Electric  field  measurements  across  the  Harang 

discontinuity,  J.  Geophys.  Res. ,  79;4620. 


Kamide  and  Akasofu  found  the  Harang  Discontinuity  to  be  the  boundary  between 
upward  and  downward  field  aligned  currents.  They  investigated  the  relationship  of 
field  aligned  currents  and  auroral  arcs  by  comparing  all  sky  pictures  with  the  mag¬ 
netometer  measurements  on  board  the  Triad  satellite  and  concluded  that:  (1)  in 
the  evening  sector,  discrete  arcs  are  in  general  associated  with  upward  aligned 
field  currents  and  (2)  no  discrete  arcs  are  associated  with  downward  field-aligned 
currents.  Figure  24  shows  the  Triad  magnetometer  traces  and  corresponding  all 
sky  picrtures  for  a  period  when  the  satellite  passed  over  the  Harang  Discontinuity, 
(shown  at  the  Station  Inuvik  on  the  map). 

3.5.4  RIOMETER  MEASUREMENTS 

A  riometer  measures  the  absorption  of  cosmic  radiation.  Riometer  measure - 
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ments  can  give  indications  of  particle  influx  prior  to  a  substorm.  Weber  et  al 
noted  the  increase  of  riometer  absorption  for  several  hours  prior  to  a  sudden  sub- 
storm  commencement  on  4  August  1972.  In  general,  during  the  austral  winter  of 
1972,  Weber  et  al's  observations  of  absorption  increases  were  associated  with 
suhstorms. 


71.  Kamide,  Y. ,  and  Akasofu,  S. -1.  (1976)  The  location  of  the  field-aligned 

currents  with  respect  to  discrete  auroral  arcs,  J.  Ceophys,  Res. , 
^:3999-4003. 

72.  Weber,  E.J.,  Mende,  S.  B. .  and  Esther,  R.  H.  (1976)  Optical  diagnostics 

August  1972  PCA  Event,  J.  Ceophys,  Res.,  81:5479-5487, 
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